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Abstract 



It is shown that there are nonlinear sigma models which are Darboux integrable and possess 
a solvable Vessiot group in addition to those whose Vessiot groups are central extensions 
of semi-simple Lie groups. They govern harmonic maps between Minkowski space M}'^ and 
certain complete, non-constant curvature 2-metrics. The solvability of the Vessiot group 
permits a reduction of the general Cauchy problem to quadrature. We treat the specific case 
of harmonic maps from Minkowski space into a non-constant curvature Lorentzian 2-metric, 
A. Despite the completeness of A we exhibit a Cauchy problem with real analytic initial 
data which blows up in finite time. We also derive a hyperbolic Weierstrass representation 
formula for all harmonic maps from M}'^ into A. 

Keywords: Pseudo-Riemanian surface, harmonic map, Cauchy problem, Weierstrass 
representation, Darboux integrability 



1. Introduction 

The linear wave equation 



on n-dimensional Euclidean space is one of the most fundamental in all of mathematics and 
physics, modeling a host of phenomena, often to a first approximation, but sometimes exactly. 
At the simplest level, n = 1, equation (1) models the motion of a frictionless elastic string 
moving in the plane under its elastic forces. In this case u{t,Xi) represents the transverse 
displacement from a reference x-axis, fixed in the plane, of the point on the string located 
at rest position Xi. It can be shown that, contrary to what is found in many textbooks on 
this elementary topic, equation (1) is not an approximate model for small oscillations of the 
string but rather an exact physical model [5], [9]. 

It is reasonable to wonder about the problem of an elastic string constrained to vibrate 
on any Riemannian surface, such as, say, the surface of a sphere, or more abstractly, on a 
2-manifold with a smooth Riemannian metric, such as the Poincare half-plane. To this end, 
let M be a Riemannian surface equipped with smooth metric g, and let P : M^^ — )■ M be a 

smooth map such that for each r, P(-, r) : R — >■ M is an immersion into M whose image is a 
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curve that models an elastic string vibrating without friction on M . Here r represents time, 
and ,^ is a parameter along the curve (string). Let V denote the Levi-Civita connection for 
(M, g). It can then be shown that the equation of motion of the string is given by [9] 

Vr.r, = c^r.r^, (2) 

where subscripts such as the 'r' in Vt denote partial differentiation with respect to r, and 
c is a constant related to the tension in the string and its line density. Equation (2) is the 
Euler-Lagrange equation corresponding to Lagrangian density 

We see that F can be viewed as a harmonic map T : R^'-*^ (M, g), where M^'^ is Minkowski 
space with metric 

In this case S{r) is the energy of F: 

If 9F°9F^ 

where we set xi = ^,X2 = r. 

Harmonic maps where the source manifold is a Minkowski space M with its usual fiat 
Lorentzian metric are generally called wave maps. (In the 2-dimensional case M = M^'^, 
they are sometimes called 1+1-wave maps.) There is a well-known geometric literature on 
wave maps that has developed over the last two decades, especially regarding their existence 
as solutions of completely integrable systems; see [13] for a textbook account with many 
references. There is also a closely related physics literature where the relevant systems are 
known as nonlinear sigma models; see [27]. 

While (1) is a linear partial differential equation and the equation of a string in the 
Euclidean plane is a linear problem, this is not the case for strings vibrating on a curved 
Riemannian — or, for that matter, pseudo-Riemannian — surface, despite the innocuous ap- 
pearance of equation (2). Indeed, if there are no elastic forces acting on the string, then the 
left-hand side of (2) vanishes and the equation is that of a geodesic in M; this reflects the 
expected result that the equilibrium conflguration of the string lies along a geodesic. As is 
well known, geodesic equations are nonlinear ordinary differential equations, and in general 
they cannot be exactly solved in terms of known functions. Thus, equation (2) precisely 
generalises the geodesic equation of a Riemannian or pseudo-Riemannian manifold, and at 
the same time, it generalises the wave equation (1). It derives its significance, in part, from 
this fact. 

Thus, just as in the case of geodesic flow, a question naturally arises as to those Rieman- 
nian or pseudo-Riemannian surfaces for which the corresponding wave map equations gov- 
erning the vibration of frictionless elastic strings exhibit "integrable behaviour." Of course, 
this question is not new: it is well known that for a wide class of manifolds M, the cor- 
responding equations are integrable in the sense of sohton theory, [13], [20]. But here we 
are interested in a stronger and more explicit form of integrability known as Darhoux in- 
tegrability [16], [21], [1], [28], [29]. Differential systems satisfying this form of integrability 
often admit Weierstrass-type representations, and, as has recently been shown, their Cauchy 
problems can be reduced to ordinary differential equations of Lie type [2], [26], [6], [7], [10]. 
In this context we note an ommision in the recent classification of Ream [17] which aims to 
give a complete list of all wave maps systems into Riemannian surfaces which are Darboux 
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integrablc at order 2 or 3. A small error in a computation caused the author to miss an 
important case which we describe herein and which, we believe, completes the classification. 

Another intriguing question in relation to harmonic maps from split signature source 
manifolds is that of global existence for the Cauchy problem. In this case, Gu [14] established 
global existence when the target manifold has a complete Riemannian metric. However, for 
split signature target metrics the global existence question appears to be more delicate. 
In that case Gu's theorem doesn't hold, as observed, for instance, in work of Terng and 
Uhlenbcck [20]. For this reason we have chosen to focus most of our attention in this paper 
on 1+1-wave maps into a certain Lorentzian 2-metric, (L, A). This metric is distinguished in 
a number of ways. We will show that the solution of the Cauchy problem can be expressed 
as the flow of a vector fleld which can be viewed as a curve in a certain solvable Lie algebra r. 
The Lie algebra r can in turn be identified as a realisation of the Vessiot algebra associated 
with the nonlinear sigma model that governs wave maps into (L, A). It follows that, because 
V is solvable, the theory of systems of Lie type permits a reduction of the Cauchy problem 
to quadrature. This makes it possible to undertake very explicit analysis of wave maps into 
this Lorentzian metric. 

The plan of our paper is as follows. In section 2 we will review the main facts about 
Darboux integrable exterior differential systems, and in section 3 we will briefiy review the 
known classification results on Darboux integrable surface metrics — that is, those Rieman- 
nian 2-metrics (M, g) such that the nonlinear sigma models for harmonic maps M}'^ — >■ (M, ^f) 
are Darboux integrable at order 2. In section 4 we will study the Cauchy problem for such 
harmonic maps into the metric (L, A); because we wish to convey a clear geometric picture 
about the nature of the harmonic maps we construct, we also include a subsection in which 
numerical solutions of an initial-boundary value problem are depicted as clastic strings vi- 
brating on the Lorentzian surface with metric {L,X). This is contrasted with the string in 
the Euclidean plane having precisely the same initial/boundary conditions. In section 5 we 
will use the theory of Darboux integrable exterior differential systems [1] and the generalised 
Goursat normal form [23], [24] to construct a hyperbolic Weierstrass representation for wave 
maps into (L, A). Finally, section 6 is an appendix containing some technical details related 
to the construction in section 5. 



2. Preliminaries 

2.1. Harmonic maps 

Let (M, g) and (A^, h) be Riemannian or pseudo-Riemannian manifolds and ip : M ^ N 
a smooth map. The energy of (p over a compact domain D C M is 

If .. d(p°'dipl^ 
4v) = 2 j^9"\x)Kp{^) ^^^ dvolM- 

The critical points of e{ip) satisfy the partial differential equation (PDE) 

where A is the Laplacian on M and A^^^ the Christoffel symbols on A^. A map </? is said to 
be harmonic if it is critical for e{ip). 

Harmonic maps generalise harmonic functions and geodesies, and have been under intense 
study since the pioneering work of Eells and Sampson [11]; sec [3] and [15] and references 
therein for comprehensive introductions. If domain((/9) = M = M, then harmonic maps are 
geodesic fiows. If codomain(99) = N = M, then harmonic maps are harmonic functions. If 
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(M, g) is pseudo-Riemannian, then harmonic maps are called wave maps or nonlinear sigma 
models. In this paper our focus is on wave map systems — specifically, the case (M, g) = 
{M.^'^,dx dy). Any further reference to wave maps in this paper means that the domain space 
M is Minkowski space R^'^ with its standard flat metric. 

2.2. Darboux integrable differential systems 

As the name implies, the notion of Darboux integrability originated in the 19th century 
and was most significantly developed by Goursat [12]. Classically, it was a method for con- 
structing the "general solution" of second order PDE in one dependent and two independent 
variables 

-^(•^) ^1 'Uyi '^xxi '^xyi ^yy) 

that generalised the so called "method of Monge." It rehes on the notions of characteristics 
and first integrals. We refer the reader to [12], [16], [21], [22] for further information on clas- 
sical Darboux integrability. There are also extensive studies of Darboux integrable systems 
relevant to the equation class under study in the works [28] and [29] . 

In this paper we use a new geometric formulation of Darboux integrable exterior differ- 
ential systems [1]. At the heart of this theory is the fundamental notion of a Vessiot group, 
which, together with systems of Lie type, are our main tools for the study of the Cauchy 
problem for wave maps. 

The PDE that govern wave maps can be put into the semilinear form 

Uxy ^ f(x,y,U,Ua:,Uy), U, f E W 

via a coordinate transformation of the form ^ — x + y,T = x — y. Each solution possesses 
a double foliation by curves called characteristics. Such PDE often model wave-like phe- 
nomena, and projection of these curves into the independent variable space describes the 
space-time history of the wave propagation. The characteristics of u^y = f are the integral 
curves of a pair of rank n + 1 distributions 

H, = {D, + Dyf ■ d^^^, , H2 = {Dy + D,f ■ d^^^, d^J , (4) 

where 

Dx = + U^du + U^xdu^ + fduy, Dy = dy + Uydu + fdu, + Uyyduy 

are the total differential operators along solutions of the PDE. Note that if 6 is the standard 
Cartan codistribution for u^y — f, then Hi® H2 — ann0. The distributions are well- 
defined with canonical structure. 

Definition 2.1. //A is a distribution on manifold M, then a function h : M ^ W is said 
to be a first integral of A if Xh = for all X E A. 

For later use, we will briefly review the key results we require from the theory of Darboux 
integrable exterior differential systems, [1]. Let Aj = anni^j, and denote by H°° and A°^ the 
final terms in the derived fiags of i^j, Aj, respectively. 

Definition 2.2. Let be a Pfaffian system on manifold M , and suppose that there exist 
auxiliary Pfaffian systems Ai, A2 such that 

a) Ai + A^ = T*M and Ai + A^ = T*M; 

b) A~ n A~ = {0}; 

c) Ai n A2 = 6. 

Then (M, 6) is said to be Darboux integrable on M, and (Ai, A2) is said to be a Darboux pair 
for 6. 
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Condition b) implies that there are no first integrals that are common to Ai and A2, 
while the crucial condition a) implies that each of Ai and A2 possess sufficient first integrals. 

For wave map equations, we have: 

Lemma 2.1. A semilinear system Uxy = f with it, / G M" is Darboux integrable at order k 
if each of its characteristic systems Hi has at least n + 1 independent first integrals order k. 

A key fact, proven in [1], is a result from the inverse problem in the theory of quotients: 

Theorem 2.2. Let (M, 6) be a Pfaffian system that admits a Darboux pair. Then there are 

Pfaffian systems {Mi,0i), {M2,02) which admit a common Lie group G of symmetries such 
that: 

1. The manifold M can be locally identified as the quotient of Mi x M2 by a diagonal 
action of G; 

2. We have the identification 



where tt, : Mi x M2 M^, i — 1,2 are the canonical projection maps; 

3. The quotient tt : Mi x M2 M to the diagonal G-action defines a surjective superpo- 
sition formula for 6. 

In this context, a superposition formula for (A/, 0) is a map tt : Mi x M2 — >■ M such that 

if ai :Ui ^ Mi are integral submanifolds of (Mj, then tt o ((ji, (T2) : Ui XIJ2 ^ M is an 
integral submanifold of 6. A superposition formula is surjective if every solution of d can 
be expressed in this form for a fixed superposition formula tt as ai range over the integral 

manifolds of di. 

Reference [1] is devoted to a proof of this theorem and to the identification and explicit 
construction of all the entities mentioned there, including the Lie group of symmetries G, 
known as the Vessiot group of the Darboux integrable exterior differential system 0. It 
is proven that the isomorphism class of the Vessiot group is an invariant of the Darboux 
integrable EDS 6. 

In this paper we shall construct the Vessiot group G and superposition formula for the 
Pfaffian system 6 governing harmonic maps R {L,X). We will also construct the 

auxilliary Pfaffian systems (Mj, di). The integration of then rehes upon the integration of 

(Mj, di). By the generahsed Goursat normal form [23], [24], we will show that for each i — 

1, 2, {Mi, Oi) is locally equivalent to a partial prolongation of the contact system on J^(]R, M?). 
In turn, this will permit us to construct a hyperbolic Weierstrass type representation for 6. 

3. Known Dcirboux integrable sigma models 

The first person to treat the Cauchy problem for wave maps into Riemannian targets was 
Gu Ghao-Hao, [14]. He established the fundamental result that for smooth initial data, wave 
maps into complete Riemannian metrics have long-time existence. Gu's work initiated many 
further investigations where regularity constraints on the initial data have been significantly 
relaxed. Furthermore, some higher dimensional problems have been treated; see [18]. 

Our interest in the present paper is in those nonlinear sigma models which are Darboux 
integrable. This question was recently investigated in the thesis of R. Ream [17] who estab- 
lished that at order 2 or 3, the only Riemannian metrics that give rise to Darboux integrable 
nonlinear sigma models are given locally by 




and g2 — 



du^ + dv"^ 



(5) 



91 = 



1 -Fe" 



1 - 
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It is obvious that these metrics can be identified by a complex translation in u, so Ream's 
is a real classification. On the other hand, it is easily shown that the metrics gi,g2 are not 
locally isometric and that the corresponding sigma models are not locally contact equivalent. 
The latter follows from the fact that the corresponding Vessiot groups are not isomorphic. 
However, it turns out that Ream's list is not quite complete, due to an error which is easy 
to repair as we show below. There is an additional metric which is Darboux integrable at 
order 2 that is locally inequivalent to both gi and g2- Before introducing it, we shall analyse 
the Ream metrics a little further. 

The signs in front of the exponential terms in (5) are isometric invariants, while the 
coefficients of u can be rescaled without changing the Euler-Langrange equations — that is, 
without changing the nonlinear sigma model. We also point out that Ream's metric gi is 
locally equivalent to the Ricci soliton metric [8] otherwise known as the "Witten black hole:" 

du'^ + dv'^ 

9^ 



Indeed, the transformation 

_r 9 _r 9 

u = e 2 cos-, v = e 2 sin - 
2' 2 

identifies gs with 

_ldr^ + d9^ 
" 4 1 + ' 

which gives rise to the same nonlinear sigma model as gi. 

It is useful to note that, whereas gi can be represented in conformal coordinates in several 

ways, it also has a representation as a surface of revolution metric obtained by taking polar 
coordinates for g^. x = rcos9, y = rsm9, and then setting r = sinhp to transform gi to 

5,4 = + (tanh^ p) d9'^. 

Thus, the nonlinear sigma model with Langangian density 

L(u) = UxUy + VxVy tanh^ u 

gives rise in spacetime coordinates (r, ^) to its equations of motion 

sinhw 

(6) 



cosh u ^ 

Vtt - = ——{UrVr - U^V^). 

cosh-u smhw 



These equations are reminicent of a system solved in [4], namely, their system (4.13). 
However, they are not quite the same. Nevertheless, the sigma model for 

du^ + (cot^ u) dv^ 

with Lagrangian density L = u^Uy + v^Vy cot^ u has, in spacetime coordinates, the equations 
of motion 

COS , o o , 

2 ^ ' 

Vtt - = ——{UrVr - U^V^). 

cosh u smh u 
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These agree precisely with equations (4.13) of [4], and in this case the authors have con- 
structed an expUcit solution in terms of arbitary functions. This system, too, is Darboux 
integrable. It can be shown that all these Darboux integrable wave map systems have Vessiot 
groups which are isomorphic to either SL{2) x R or to SO{3) x R. 



4. A nonlineeir sigma model with a solvable Vessiot group 

4.I. Completing the Ream classification 
The surface of revolution metric 



du'^ + ^dv^ (8) 



arose in [25] as an application of hyperbolic reduction. It was shown that wave maps into 
(8) can be constructed as integrable extensions of the scalar hyperbohc PDE 



u 



xy 



U 



l-ulJl-ul (9) 



and that the Vessiot group of the integrable extension is 4-dimensional and solvable. Equa- 
tion (9) occurs in Goursat's well known classification of Darboux integrable scalar second 
order PDE in the plane. Moreover, it is straightforward to show that the nonlinear sigma 
model for (8), 

Uxy H f = 0, V^y - ^ = 0, (10) 

is itself Darboux integrable [25]. 

Lemma 4.1. The metrics 

dx^ + dy"^ dx'^ + dy"^ ,2 1,2 

9i = , , _^ , 92 = — gp = dx + —dy 



1 + ' ^ 1 - 



X 



are pairwise inequivalent. (The metrics gi,g2 can, however, be identified by the complex 
transformation x ^ x + 7r-\/^. j 

Proof. This is a direct application of the Cartan equivalence method for 2-metrics. Firstly, 
the metric gp — dx^ + dy^/x^ has orthonormal coframe cu^ — dx, lo^ — dy/x and lifted 
coframe 

© = {6*^ = cos Au;^ - sin Au;^, 6*^ = sin Au;^ -\- cos Au;^}. 
The connection form on the principal 5'0(2)-bundle is 

= dA + -w^, 
x 

giving rise to the canonical structure equations 

de^^-u^^^9\ 
de^^uj^^A9\ 

x^ 



The coframe {6^, 6'^, cu^^} has rank 2 and order 1, and its invariants are spanned by the 
functions K and Ki, where K denotes the Gauss curvature of the metric, and the functions 
Ki, K2 are defined by the equation dK = Ki6^ + K2 9^. For the metric gp, v/e have 

2 4 cos A 4 sin A 

= 2^ -^1 = ^> -^2 = , 

x"* x"* 



and these functions satisfy the relation 



(11) 



Performing the equivalent calculation for metric gi of (5) gives 

3 



and for g2 of (5) we have 



1 e 2 

(12) 



fi 2 

(13) 



2 + 4 

Since these functional relations are different, it follows from Cartan's solution of the local 
equivalence problem that the metrics are pairwise inequivalent. □ 

Corollary 4.2. The nonlinear sigma model (10) for wave maps into the metric 

gp — du^ H — ^rif 2 

extends the Ream classification [17] of metrics for which the wave map equation is Darhoux 

integrable at order 2 (or 3). 

Remark 4.1. The author in [17] missed the metric gp in his classification because of an 
allowable form of an integating factor with a free parameter C in which he inadvertently 
ommitted the value C = in the prelude to equation (4-55) on p53. Setting C = prior 
to (4-55), following his subsequent calculations and then making an elementary change of 
variable leads to metric gp. 

4-2. The Cauchy problem 

According to Gu's theorem [14], the Cauchy problem for wave maps into complete Rie- 
mannian metrics has global existence for smooth initial data. The situation appears to be 
rather different for wave maps into split signature metrics. For instance, Terng and Uhlen- 
beck discovered a bi-invariant Lorentzian metric on SL{2,M^) and showed that there were 
smooth, finite energy wave maps which blew up in finite time. 

With the ultimate aim of studying the properties of harmonic maps between pseudo- 
Riemannian manifolds, we will in the remainder of this paper turn our attention to wave 
maps into Lorentzian metrics. Metric (8) can be put into the conformal form 

du'^ + dv"^ , , , 
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via the transformation 

u = \/2^, V = V. (15) 

In this section we will show that the general Cauchy problem for all harmonic maps from 
Minkowski space M^'^ into the semi-Riemannian metric 

duj - dul 

is reducible to quadrature and that, given any initial data, its solution can be constructed 
by solving the initial value problem for an ordinary differential equation of Lie type [6], [7], 
[10], [26] associated to the free and transitive action of a solvable Lie group on M*^. 

The nonlinear sigma model in this case is also Darboux integrable. In lightcone coordi- 
nates on both source (x, y) and target (m, v) spaces 

{u^,U2) ^ ( — 2 — ' — 2 — ) " ^^'"^^ 

it has the form 

Each of the characteristic systems i?i and in (4) have 4 independent first integrals, 
which we shall label 

y, ^1, ^2, /^s for ifi and x, cti, q;2, as for B.^- 

Specifically, we have 

_ UyVy n _dPl Vyy Uy 

Pi — \ ) P2 — P3 — , , 

u + V ay Vy u + V 

and 

UxVx dai Uxx Vx 

Oil = , , Oi2 = — , q;3 = ■ — . 

u-\-v ax Ux u -\- V 

Adapting these first integrals as new coordinates on the differential equation submanifold of 
J^(]R^,]R ) considerably simphfies the Pfaffian system for (17). We obtain: 

= { dbi — 62%, dai — a2dx, dzi — z^dx ^{zi + Z2)dy, 

dz2 ^{zi + Z2)dx — z^dy, dz^ — (oi + zza^)dx -—dy, 

Z3 Z4 



dz4 ^^dx — (63^4 + bi)dy }, 

Z3 



(18) 



where 



zi = u, Z2= V, Z3 = Ux, Z4 = Vy, ai = ai, bj = /3j. 
Our main results are contained in the following theorem: 
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Theorem 4.3. Consider the initial value problem for the harmonic map system 



u+v u+v 

with initial data 

du dv 

u\,=4>,. .,, = 02, ^,^=^1, a^|^=^2' (19) 

where 'j is a curve with tangents nowhere parallel to the x- or y-axes, n is a unit normal 
vector field along 7 and (f)i,ipi are smooth functions along 7. 

1. The Cauchy problem has a unique smooth local solution. Moreover, the unique local 
solution is expressible as the solution of an ordinary differential equation of Lie type 
associated to a local action of a solvable Lie group on M . 

2. Given the unique local solution {u,v) from part 1., the Cauchy problem for harmonic 
maps 

(M^'S dxdy) [M, 



du\ — du\ 
is given by 

Ui—U-\-V, U2 — U — V, 

where Ui , U2 satisfy initial conditions 



dui du2 
— = -01 + -02, ^ 
on on 



^^11^=01 + 02, ^i2|T = 01-02, ='01+'02, =V'l-V'2- 



Proof. By hyperbolicity, the problem is locally well posed. Let ki{y),k2{y) be arbitrary, 
smooth real- valued functions, and consider the overdetermined PDE system defined by (17) 
together with the additional PDE 

Pi = h{y), (32 = k{y), /33 = k2{y), (20) 

where the dot denotes y-difi^erentiation. It can be shown that this overdetermined system S' 
is involutive and, moreover, admits a 1-dimensional Cauchy characteristic distribution. Now 
suppose we fix a smooth curve 7 embedded in a portion of the xy-plane J\f, and suppose that 
Cauchy data is prescribed along 7 as described in (19). Then by an argument similar to ([21], 
Proposition 3.3), 7 can be lifted to a unique curve 7 : / C M ^ J^(A/', M^) which agrees 
with the Cauchy data. Let i : 'HkiM ~^ J'^iM,^^) denote the submanifold in J^(A/', M^) 
defined by PDE system £' ^ and let denote the contact system on J^. Let Oi^^^^^ = l*Q 
be the Pfaffian system whose integral submanifolds are the solutions of S'. We now claim 
that the functions ki, k2 can be uniquely chosen in such a way that 7 can be extended to a 
1-dimensional integral submanifold 7 of dki,k2- To see this, consider the Pfaffian system 

9ki,k2 — |<^^ — du — Uxdx — Uydy, cu'^ = dv — v^dx — Vydy, 

oj^ — dux — Uxxdx ^^dy, 00'^ — duy ^— ^o?x — Uyydy, (21) 

u+v u+v 

oj^ — dvx — Vxxdx ^^dy, oj^ — dvy ^^—^dx — Vyydy\. 

u + v u + v "«J 

(22) 



10 



Pulling this system back by 7, wc observe that a;-'", u;^ pull back to zero by construction. The 
forms cu'^ and cj^ define the functions Uyy and Vyy along 7, while 00'^, 00^ define the functions 
Uxx and Vxx along 7. All these functions are expressed in terms of the Cauchy data 4>i,ipi. 
Substituting these expression back into (20) uniquely determines the functions ki, k2 in terms 

of (f>i,iJi. 

We will now use the first integrals /3j recorded above to demonstrate that ,^2 has 
a one-dimensional Cauchy characteristic distribution and that, in particular, the Cauchy 
characteristic vector field can be chosen to be a curve in a certain Lie algebra — the Vessiot 
algebra [1] of system (17). It will be seen that the Cauchy characteristic vector field is gener- 
ically transverse to the Cauchy data and extends the one-dimensional integral submanifold 
of 0ki,k2 to the solution of the Cauchy problem. Because the Cauchy characteristic vector 
field is a curve in a Lie algebra, this extension from a one-dimensional to a two-dimensional 
integral of 0ki,k2 is an ordinary differential equation £, of Lie type. Its coefficients and initial 
conditions are fixed by all the data present in the problem, including the Cauchy data. Any 
solution of £ (independently of its initial conditions) permits a Lie reduction of £ and will 
permit us to solve the I VP for £,. 

We find by explicit calculation that the vector field 



Ck,M = dy-Ri- k2{y)R2 + /ci(|/)i?3 (23) 

spans the Cauchy characteristic vectors of dk^,k2: where r = i?2, -R3, -R4} is a basis for a 
solvable Lie algebra with nonzero structure 

[Ri, R2\ — Ri, [Ri, Rs] — —Ra, [R2, -R3] = Rz 

with local expression 

R^ = -Zid^^, R2 = -Z4d^^, R3 = ^^-t^^^i + —9^3 + 9^4, i?4 = - 9^2. (24) 

Z4 Z4 

Note that r generates a free transitive local Lie group action on M^. 

We have shown that, given Cauchy data along a generic smooth curve 7 in the 

xy-plane, there is a unique choice of ki and ^2 such that the Cauchy data can be lifted to 
a 1-dimensional integral submanifold of OkiM- Indeed, following the outlined procedure in 
the case of the important curve 7(x) — (x, x) gives 

, I 2V>i-02 + (klxhx - V^{<Plxi^2 + -0102a;) 

4(0, + 02) 

and 

(01 + 02)(202.. - 2^V2x) + (\/2(^l - V'2) - (01x + 02.))(02. " V^V'2) 



k2{x) = 



V2(0l + 02)(V202a:-2V'2) 



Thus for any smooth, real- valued functions ki{y),k2{y), Ok^^k2 admits a Cauchy charac- 
teristic vector field ikiM- Furthermore, we have shown that the Cauchy data along 7 can 
always be lifted to a unique 1-dimensional integral manifold of 0^1, fc2 for a unique choice of 
/ci,/c2, depending upon the Cauchy data. The Cauchy characteristic vector field can then 
be used to extend this 1-dimensional integral manifold to a 2-dimensional integral mani- 
fold, provided that the Cauchy characteristic vector field is everywhere transverse to the 
1-dimensional manifold. But because we have chosen initial data which has tangents par- 
allel to neither coordinate axis in the Minkowski plane, transversality is guaranteed; hence, 
the Cauchy characteristic vector field ^kiM provides an extension of the initial data integral 
manifold 7 to the unique smooth local solution of the Cauchy problem. 

To see that the extension by S,ki,k2 from a 1- to a 2-dimensional integral manifold is an 
ODE of Lie type, it is sufficient to note that ikiM is a curve in the Lie algebra r. □ 
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Remark 4.2. At this point, the occurrence of Lie algebra v may appear somewhat mysterious. 
Its geometric origin and construction is discussed in the Appendix in section 6 and used again 

in section 5. 

Lemma 4.4. The metric 

du^ — dun 

^ = — 

2ui 

is complete. 

Proof. In lightcone coordinates {u,v) on {L,\), the geodesic equations are 

u — , V 



u-\-v u-\-v 
The complete set of geodesies are 



u 



^ (^^{ai + 6i)e(""''^* + ah{h - a)(ai + hi)t - a% - 2abai + aib'^^ , 

V = j^^^^-^ (^b\ai + 6i)e-("-'')* - ab{b - a){ai + bi)t + a% - 2abbi - aife^^ , 

where a,b,ai,bi are constants depending upon initial conditions. Note that u = constant 
and V — constant are geodesies of A. □ 

Example 4.1. We will use the above solution of the Cauchy problem to show that, even 
with real analytic initial data, harmonic maps from Minkowski space into the Lorentzian 
metric A may develop shocks in finite time, despite the fact that the metric is complete, in 
contrast to Gu's theorem [14] . 

Consider the initial value problem for (17) subject to the initial conditions 

On Ov 

»,, = 1--. ",,=2-. snu = 0, ^^^ = 0, (25) 

where 7 is the curve y = x in the xy-p\ane and n = 2~^{dx — dy). We lift 7 to a 1-dimensional 
integral of 6ki,k2 ^^r suitable unique functions ki,k2, as follows: observe that 



Hence we have the partial hft 



(^1,^2,^3,^4)1^= y--x, 2x, 1 j . 

We can extend this partial 1-dimensional integral by direct substitution, and we find that 
the curve cr : M — >• HkiM defined by 

/ 112 

{zi, Z2, Z3, Z4, ai, a2, 03, bi, 62,&3)|^ ^ [l - x,2x,--, 1, 



2' ' 2 + 2x' (2 + 2x)2' 
1 2 1 



2 + 2a;' 2 + 2a;' (2 + 2a;)2' 2 + 2a; 
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along y — X is a. 1-dimensional integral manifold of OkiM where 
Hence the Cauchy characteristic vector field of 6ki,k2 is 

Computing the flow for (26) subject to the initial conditions above along 7 gives the explicit 
solution of Cauchy problem (25): 

{-yy/l + x - y/TTy + xy/TTy - V^ + x + (1 + yf^"^ 
4V1 + 2/ - Wl + x + y^Jl +y - AyVl + x + xyTTy 



V — 



The corresponding solution ui — u+v, U2 — u—v for the metric (16) in space-time coordinates 

^ = y = 



IS 



(5C - 5t + 10) V4 + 2^ + 2r - (8 + 4^ V4 + 2^ - 2t 
~ 2V4 + 2^ - 2t 

(3t - 3^ - 6)V4 + 2e + 2T + 8V4 + 2^ - 2t 



U2 = 



2V4 + 2^ - 2t 
That is, the functions Mi(^,T) satisfy 



dr^ de 2ui\\dT ) ^ \ dT ) \ J \ 

d'^U2 d^U2 1 / dui du2 dui du2 



i2 



(27) 

Note that in space-time coordinates, the initial conditions have the form 

M^,o) = i + l^, «2(C,o) = i-^e, ^(^'0) = ^(^'0) = 0- 

Despite the analyticity of this initial data, the solution exhibits a shock wave whose image 
is moving in the ^ • r-plane along the curve 
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At a given value of { > 0, both solutions ui, U2 experience infinite discontinuities at a time r 
soon after t = ^ + 3^^. This behaviour is excluded in the case of harmonic maps from 
Minkowski space to complete Riemannian metrics according to the result of Gu [14] . 

Example 4.2. On the other hand, if we take a slight variation of these initial conditions 
and choose 

ui^ = l- x"^, VI 

we find that 

hiy) 

Then computations similar to those of the previous example show that the Cauchy charac- 
teristic vector field 

dy-Ri^^R2- y'^Rz (29) 
of 6_y2^i/y leads to the simple globally defined solutions: 

^^l = l-^^V^ U2 = l-l{e + r') (30) 
satisfying the following initial conditions in space-time coordinates: 

^l(e,0) = l, U2{t0)^l-le: ^(e,0) = ^(e,0)=0. 

Note that in both Examples 4.1 and 4.2, the energy of the initial data is infinite. However, 
in the latter case the energy density is an analytic function of the parameter ^, whereas in 
the former it has a simple pole in ^. 

4-3. The Cauchy problem and systems of Lie type 

One of the reasons we are able to get compact, simple formula for solutions of the Cauchy 
problem is that the Vessiot group G (whose Vessiot algebra is r) is solvable. As we have 
shown, the solution of the Cauchy problem is equivalent to the flow of a specially constructed 
vector field ^kiM that depends on the choice of initial conditions (f)i,ipi and their 2-jets. We 
have observed that C,ki,k2 is a curve in the Vessiot algebra r of the nonlinear sigma model 
(17). That is, the differential system that explicitly determines the solution of the general 
Cauchy problem is an ODE of Lie type. It is well known [10] , [6] , [26] that if, as in this case, 
the group action associated to the system of Lie type is solvable, then the solution can be 
reduced to quadrature. Thus, we have 

Proposition 4.5. The solution of the Cauchy problem for wave maps of Theorem 4-3 

(M^'\ dxdy)^(M,^^^^=^) 

is reducible to quadrature. The solution can constructed by extending the Cauchy data to a 
1-dimensional integral manifold of 0ki,k2 then flowing this out to the unique solution via 
the flow of the Cauchy vector ^kiM (of Ok^.k-z)} where ki{y),k2{y) are real-valued functions 
that are uniquely determined by the initial data. 



du 



dv 



7 = ^ =0, — =0, 



(28) 
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4-4 ■ Numerical example 

Recall from the Introduction that an clastic string vibrating without friction and con- 
strained on a surface equipped with a Riemannian or pseudo-Riemannian metric g is math- 
ematically modelled by wave maps into g. To gain additional insight into the nonlinear 
dynamical systems under study, we include a brief section on the interpretation of nonlinear 
sigma models as partial differential equations that precisely describe this motion. We are 
interested in particular in the nonlinear sigma model for harmonic maps from Minkowski 
space M}'^ into the metric A. In order that boundary conditions also be included, we have 
computed and graphed numerical solutions of two wave map systems. Figures 1 and 2 each 
display graphs of wave maps subject to the boundary and initial conditions 



ui{0,r) = 0, Ki(2,t) = 2, K2(0,t) = 1, U2{2,r) = 1, 




0.5 1 1.5 2 0.5 1 l.S 2 




Figure 1: Evolution of the string in Euclidean plane. 




0.5 1 1.5 2 0.5 1 1.5 2 



Figure 2: Evolution of the string in Lorentzian metric A. 

The values of U2 are plotted along the vertical axis and those of ui along the horizontal 
axis. Initially, the string occupies < -Ui < 2, ^2 = 1; is clamped at {ui,U2) = (0,1) and 
{ui,U2) = (2, 1), and the ?i2-componcnt has initial velocity Ui{ui — 2) at Ui. Figure 2 plots 
the time (r) evolution of curves u : M^^"^]. — )■ {L,X) for r G [0,T], for some T > 0. For 
comparison, exactly the same boundary and initial conditions are imposed on wave maps 
into the Euclidean plane and displayed in Figure 1. Note that 1x2 = 1 is not a geodesic of A, 
so even without the initial downward initial velocity the string would vibrate from its initial 
position. However, the nonzero initial velocity increases the amplitude of the vibration and 
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leads to more interesting behaviour. By the fourth frame in Figure 2, the curve/string ceases 
to be immersed at Ui = 1 and at some epochs during its motion. It also loops around to 
intersect itself briefly. At this time, symmetric "waves" form on each side which persist to 
the end of the computation. The lifetime of this solution is not known to the authors.^ 

5. Hyperbolic Weierstrass representation 

In this section we use the Darboux integrability of the wave map equation (17) to compute 
its general solution and hence construct a hyperbolic Weierstrass-type representation for 
wave maps into the corresponding metric. According to [1], we pull back to suitable 

integral manifolds Mi,M2 of h[°°^ and i72°°^ respectively. It is convenient to define Mi by 
y — bi — b2 — bs — and M2 by x — ai — a2 — — 0. This gives Pfaffian systems 

61 — < dai — a2dx, dzi — z^dx, dz2 ^(-^i + Z2)dx, dz^ — (oi + a^Zz)dx, 

I 

aiZ/^ 
dz4 dx 

and 

02 = s dbi — b2dy, dzi — —{zi + Z2)dy, dz2 — z^dy, dzs-^^^^dy, 
I ^4 ^4 

dz4, - {bi + b3Z4)dy \, 
each of rank 5, on 8-manifolds (Mj, di). Locally Mi has coordinates 

X, Zi, Zi, Z2, Zs, Z4, Oi, 02, Qs, 

while M2 has local coordinates 

y, zi, zi, Z2, Z3, Z4, bi, 62, Ca- 
using these local formulas, we define a local product structure 

(^Mi X M2,^ie^2) , 

where Mi = Mj and locally Mi has coordinates 

X,qi,qi,q2,q3,q4,ai,a2,a3, 

while M2 has coordinates 

X,Pl,Pl,P2,P3,P4,bl,b2,b3. 

The Pfaffian systems 61 and 62 are locally generated as 



^See [9] for further discussion and examples of string motion on Riemannian surfaces. 
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and 



01 = dtti — a2dx, dqi — q^dx, dq2 -{qi + q2)dx, dq^—^ai + azq^)dx, 

, am 

dq4 dx 

Q3 



02= { dbi - 62%, dpi - — (pi + P2)dy, dp2 - pdy.dpz - ^-^dy, 

Pa Pa 



dp A - {hi + b3P4)dy 

As described in [1], the relationship between 6, 61 and 62 is that every integral manifold 

of is a superposition of integral manifolds of Oi and 62. The superposition formula is the 
map^ 

TT : Ml X M2 ^ M 

defined by 

n{x,q,a;y,p,b) = 

-qi + qiPs + QaPi - g2 + ^2^3 + QiQa qAPi + qAP2 + g2P3 - Pi - P2 + P3P2 
54' P3 ' 

?3(?4 - 1 + Ps) P4(g4 - 1 + P3 



?4 P3 



, X, 7/, a, 6 



(^^1, 2;2, ^3, ^4, X, y, a, b^. 



The usefulness of this factorisation of the integration problem for 6 is not only that the 
integration of 6^ relies on ODE while that of relies on PDE, but also that the 0i are locally 
equivalent to prolongations of the contact system on J^(]R, M^). To see this, we turn to the 
characterisation of partial prolongations of such contact systems provided by [23], [24], which 
also provide simple procedures for finding the equivalence. To implement this, we compute 
the annihilators 



ann Oi — Hi — 

\ dx + qzdqi + —(?! + 52)^92 + («3g3 + 01)^93 + ^^^94 + «25ai, 5a2, ^as f , 
L ?3 ?3 



^See section 6 (appendix) for further details on the superposition formula; in particular, how it is defined 
and constructed. 
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aim 02 = H2 



\dy + — (pi +P2)dp^ +P4dp2 + — + {hP4 + bi)dp^ + b2db^,db2, dbA . 
K P4 Pa J 

The construction we have just described can be summarised by the following commutative 
diagram, where tj, i2 are inclusions and 7ri,7r2 are projections such that tt^^j = Oj- 

(Ml X M2,^i ©^2) 





(M2,02) 



(M,0) 

Locally, the maps tti : Mi — >■ Mi and 7r2 : M2 — >■ M2 are given by 

7ri(p,a,a;) = {z,a,x), 7r2{p,b,y) = {z,b,y). 

We now show that each of the Hi is locally equivalent to the partial prolongation C(0, 1, 1) 
of the contact distribution on J^(]R, M^) — that is, the contact distribution on J^(]R, M^) par- 
tially prolonged so that one dependent variable has order 2 and the other order 3, with 
canonical local normal form 

c{o,i,i) = [dt + zld,i + zld,i + zld,2 + zld,. + zld,2, 

Let I> be a smooth distribution on the manifold M, and assume that T> is totally regular in 
the sense that V, all its derived bundles, and all their corresponding Cauchy bundles have 

constant rank on M. Denote = dimP^'^ = dim Char and xl-i — dim Char T^f^^, 
where 

Char P^i = ^^'"^^ n Char V^^. 

Below k denotes the derived length of T>. 

According to Theorem 4.1 of [23], a totally regular distribution T> on smooth manifold 
M is locally equivalent to a partial prolongation of the contact distribution on J-'^(M, R*) for 
some q if and only if 

1. The integers mi,x' ,Xi-i satisfy the numerical constraints 



2mj — m,j^i — 1, 0<j<A; — 1 



x\-i = rrii-l, l<i<k-l 
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2. If ruk — rrik-i > 1, then a certain canonically associated bundle called the resolvent is 
integrable.^ 

A pair {M.I)) that satisfies these conditions is said to be a Goursat manifold or Goursat 
bundle. Moreover, if P is a Goursat bundle on M, then it is locally equivalent to a partial 
prolongation with ~ Xj-i dependent variables at order j < k and — m^-i dependent 
variables at highest order k. This uniquely identifies the partial prolongation associated to 
a given Goursat manifold. Before discussing the question of constructing equivalences, let us 

solve the recognition problem for the distributions H^. We will demonstrate the procedure 
for Hi, which we denote temporarily by K. 
We find that 

Char K = {0}, Char K^'^ = {da,, da,}, Char i^f^ = {da„ da„ da„ d^,}. 

Char = {da,, da,, da„ dg^, d^}. 
Calculation shows that the dimensions of the derived bundles are 

dim ^ = 3, dim K^^'^ = 5, dim K^'^^ = 7, dim K^'^ =8, s > 3. 

Hence the derived length is k — 3. Below we check the first condition of a Goursat bundle: 



j 


mj 


mj_i — 1 


2mj — rrij^i — 1 




X' 





3 




6-5-1=0 







1 


5 


2 


10 - 7 - 1 = 2 


2 


2 


2 


7 


4 


14 - 8 - 1 = 5 


4 


5 



Table: Checking the numerical constraints satisfied by {M,V). 

Therefore, i^' is a Goursat bundle with k = 3, rUk — ruk-i = 1, and the only nonzero 
difference X'' ~ x]-i being at order j = 2: x^ — x^ = 5 — 4 = 1. Hence there is one variable 
of order 2 and one variable of order 3. This solves the recognition problem, and we can 

assert that K := Hi is locally equivalent to C(0, 1, 1). Next, we show how to construct an 
equivalence. Given a Goursat bundle, an efficient method for constructing an equivalence 
map was worked out in [24]; it relies on the filtration induced on the cotangent bundle. 

Denote by S*^-^^ and the annihilators of Char K^^^ and Char respectively. Then 

we obtain the filtration 

5(2) ^ 5(2) ^ 5(1)^ 

expressed explicitly as 

{dqi,dq2,dq4} C {dqi,dq2,dq4,dx} C {dqi,dq2,dq3,dx,dq4,dai}. 

The construction proceeds by building appropriate differential operators and functions. Be- 
cause rUk — ruk-i = 1, condition 2 of the definition of Goursat manifold is vacuous. Instead 



^In the original formulation of Theorem 4.1 in [23], the integrability of Char Vy}^ is an additional 
hypothesis to be checked. This is a simple task, but unnecessary, since it is easy to see that this bundle is 
always integrable and that this hypothesis can be omitted. 
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we fix any first integral of Char V^'^ denoted t, and any section Z ofT> such that Zt — 1. 
Then, define distributions n*^ inductively as follows: 

n^+^ = [z, n^], = char vl, i<e<k-i. 

There is a function ip'^ which is a first integral of II'^ such that d(p'' A dt ^ 0. The function 
(p'' is said to be a fundamental function of order k. The space of fundamental functions of 
lower order can be constructed from the filtration above by taking quotients. Specifically, as 
noted above, in this case the only fundamental functions of less than maximal order 3 are of 
order 2. They are described by the quotient bundle 

Ef^/E('^ = {dx}. 

Without loss of generality, we can take (p'^ — x to he a fundamental function of order 2. The 
construction of the equivalence map is now as follows. The function t is the "independent 
variable," and successive differentiation gives the higher order variables 

We now implement this. The first integrals of Char X^^^ are spanned by q'l, ^2, 93, and any 
function of these can be chosen to be t. If we choose (say) t — qi, then for Z we choose 

Z= —X, 

Q3 

where X is the first vector field in the basis for K above; for then Zt = 1, as required. We 
then construct the integrable distribution 11^ as described above and discover that its first 
integrals are spanned by 

qi + q2 

Hence the fundamental function of highest order (three) is 

qi + q2 

The data 

94 



t = qi, Co =x, Co = 

91 + 92 

and differentiation by Z now construct the local equivalence ip identifying K — Hi and 
C(0, 1, 1). The local inverse cxi = ■j/'f^ • ~^ -^i determines the integral submanifolds of Hi. 

An exactly analogous calculation holds for H2, and one arrives thereby at an explicit 
map o"2 = "^2^^ • ^2 representing the integral manifolds of H2- The explicit integral 

manifolds of 6 are a superposition of those of Hi and H2'- 



X 



^ 77(^^1-' (K),V'2"' 



The local diffeomorphisms ipi are guaranteed to have a block-triangular structure and are 
therefore usually not difficult to invert in practice. 

In this way we obtain remarkably compact representations for wave maps into this metric: 
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Theorem 5.1 (hyperbolic Weierstrass representation). For each collection of twice con- 
tinuously dijferentiable real valued functions k{s), h{s), m{t), fit) of parameters s,t, the 
functions 

X = f(t), y^h{s), 

mks + mk + k + k^ — tmk 



u 



mk 

+ tkrh + mk + m — ksm 



km 

define harmonic maps 

{R^''^, dxdy) ^ In, — ^- -\ by ui^u + v, U2^u-v. 

6. Appendix: The superposition formula 

For completeness, in this appendix we make a remark on the construction of the super- 
position formula tt : M\ x M2 — > M. A general construction valid for any decomposable 

exterior differential system was worked out in [1]. Below we present results for the wave map 
system studied in this paper. The hyperbolic structure of the wave map system in adapted 
coordinates is given by 

where 

Hi= \d,, + z^dz-, + — (^1 + Z2)dz^ + (a3^3 + ai)dz^ + ^^d^^ + a2da^,da2,da3 r , 
I Z3 Z3 

H2 = \dy + —{zi + Z2)dz^ + Zid^^ + —9^3 + (632:4 + &i)9^4 + b2db^,db2, 
L Z4 Z4 

Calculation shows that the infinitesimal symmetries of Hi which are tangent to the level sets 
of all the first integrals x,y,a,b, known as the tangential characteristic symmetries of Hi, 
are spanned by 

ei = i -Z4dz^,-Z4dz^, ^^-^5^1 + —5^3 + d^^, d^^ - d^^ 

y Z4 Z4 

Similarly, the tangential characteristic symmetries of H2 are spanned by 

t2 = \ zsd^-^zsd^^, ^^-^d^^ + 9^3 + —9^4, -d^-, + d^^ 

L Z3 Z3 

The structure of these Lie algebras are in reciprocal relation; namely, for Ci the nonzero 
Lie brackets are 

[61,62] = -ei, [61,63] ^ 64, [62,63] = -63, 

while for 62 we have 

[/i)/2] = fi, [/ij/a] = — /4, 1/2, /s] = fs, 
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and 

[ei,fj] =0, V ij. 
Recall the local product manifold Mi x M2: 

Ml X M2 = {q,a,x]p,b,y) 

where 

p = Pi,P2,P3,P4; q = qi, (12, q?>, (lA- 

The diagonal action : G x Mi x M2 ^ Mi x M2 takes the form 

fj.D{g){mi,m2) = (^fj.i{g'^) {nil); fj.2{g) {1712)^ , ^ g e G, (mi,m2) G Mi x M2, 

where is the action generated by Tf*ti on Mj. This is the diagonal action mentioned in 
Theorem 2.2. 

The Lie algebras Ci, C2 define the isomorphism class of the Vessiot algebra for this Darboux 
integrable system, in the terminology of [1] . Each Cj is isomorphic to a solvable Lie algebra 
Q, and each frames a neighbourhood of a point of M^. By the converse of Lie's second 
fundamental theorem, there is a local Lie group (G, m), where m : G x G ^ G denotes 
group composition, such that coincide with the infinitesimal left and right translations on 
G. From the expressions for either Ci or t2, by computing flows or otherwise, we can compute 
the function m. The components of m coincide precisely with the first 4 components of 
the superposition map cr, where we interpret {{pi}f^i, {qi}f^i) as local coordinates around 
(e, e) E G X G; e being the identity in G. Thus, we see that the Cauchy characteristic vector 
field constructed in Theorem 4.3 is a curve in the Lie algebra of tangential characteristic 
symmetries of Hi. There is an entirely analogous construction in which the first integrals aj 
of H2 are used instead of those of Hi, and in that case the Cauchy vector can be expressed 
as a curve in the tangential characteristic symmetries of H2. 
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